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双 参 数 Rayleigh 分 布 环境 因子 的 经 验 Bayes 估计 及 应 用 * 


胡 俊 梅 ， 师 义 民 ， 高 WE, EDR 
(西北 工业 大 学 应 用 数学 系 ， 西 安 710072) 


摘 要 : 本 文 针对 Rayleigh 分 布 位 置 参数 已 知 的 情形 ， 给 出 了 Rayleigh 分 布 环境 因子 的 极 大 似 然 估计 和 
经 验 Bayes 估 计 ， 并 将 环境 因子 的 估计 结果 应 用 于 Rayleigh 部 件 的 可 靠 性 评估 ， 给 出 了 该 部 件 
可 靠 度 函 数 与 失效 率 的 估计 。 最 后 的 随机 模拟 例子 表明 ， 经 验 Bayes 估计 优 于 极 大 似 然 估 计 ， 并 
且 在 考虑 环境 因子 的 情形 下 ，Rayleigh 部 件 可 靠 性 指标 的 估计 优 于 未 考虑 环境 因子 时 的 估计 。 
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在 可 靠 性 寿命 试验 中 ，Rayleigh 分 布 由 于 其 良好 的 单调 失效 率 性 质 ， 受 到 了 许多 科技 工作 
者 的 极 大 关注 ， 它 被 广泛 地 应 用 到 医学 、 通 讯 工 程 、 机 械 设计 和 航空 航天 等 领域 中 四 。 双 参 
数 Rayleigh 分 布 的 分 布 函数 和 概率 密度 函数 分 别 为 


(2-4)? 
F(r;u,0)-—1-—e^ 207^ 3 (1) 
d — LG 
g? 


f(z; p, 0) = : (2) 


式 中 z p0, 1 为 位 置 参数 ;co > 0 为 尺度 参数 。 

已 有 一 些 文献 对 Rayleigh 分布 做 了 讨论 。 文 献 [L23] 讨 论 了 单 参数 Rayleigh 分 布 的 参数 
估计 问题 ， 文 献 [3] 讨论 了 双边 定数 截 尾 情形 下 单 参数 Rayleigh 分布 的 Bayes 估 计 和 预测 问 
题 ， 文 献 | 各 讨论 了 双 参 数 Rayleigh 分 布 参数 的 条 件 置信 限 ， 文 献 [5,6] 分 别 讨论 了 Rayleigh 3€ 
fil Rayleigh Ñi -o Am, XT Rayleigh 分布 环境 因子 的 讨论 尚未 见 到 。 在 工程 实际 中 ， 环 境 
因子 作为 一 种 折合 因子 ， 可 以 扩大 样本 量 ， 提 高 可 靠 性 评估 水 平 。 为 此 ， 本 文 讨论 双 参 
数 Rayleigh 分 布 环境 因子 的 估计 及 在 可 靠 性 评估 中 的 应 用 。 


2 Rayleigh 分 布 环境 因子 的 定义 


关于 环境 因子 的 讨论 建立 在 以 下 假设 基础 上 [7]。 
假定 A 在 不 同 环境 下 产品 的 寿命 分 布 属 于 同一 Rayleigh 分 布 ， 只 是 在 分 布 参数 上 存在 差 
异 。 即 环境 I 和 环境 开 下 产品 的 寿命 分 布 分 别 为 


CET 
F(z; mo) =1~e 2 , i-1,2 
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假定 B 在 不 同 环境 下 产品 的 失效 机 理 不 变 ， 反 映 在 Rayleigh 分 布 中 ， 位 置 参 数 不 变 。 即 
pi — ua — p. 


假定 C 产品 的 残存 寿命 仅 依赖 于 已 积累 的 失效 和 当前 的 环境 应 力 ， 而 与 累计 方式 无 关 ， 
即 
Fi(z1; 4,01) = Fa(z2; 4,02) €» (21 — 1)2/201? = (£2 — H)’ /2027. 


基于 以 上 三 个 假定 ， 可 以 导出 Rayleigh 分 布 的 环境 因子 为 
k = 02/01. (3) 


FÜSET AER SE uo EXEC E. 


3 ”环境 因子 的 极 大 似 然 估计 


4 已 知 时 ， 讨 论 无 替换 定数 截 尾 模型 下 环境 因子 大 的 估计 。 

B ta < tüz L-e < Zlimi) NK H Rayleigh 分 布 容 量 为 ni 的 随机 样本 中 前 rr 个 最 小 
观察 值 (为 方便 起 见 ， 可 将 zth 的 下 标 数字 省 略 括号 ， 下 文 的 xij 表示 第 ;个 最 小 观察 值 )。 
$ ni = (zii, Tiz, t, Cim): 由 文 [8] Al, 基于 样本 r; 的 似 然 函 数 为 


1 Ti - 
Li(zi|o;) = eos lI f» e) | E= (siu, 9i). (4) 
1 2/* j=1 
将 (1), (2) 式 代 入 (4) 式 ， 可 得 
Li(zj | ci) = Aie— Qi/20 Jg 2m, (5) 


其 中 


nil [a (zij; — p) 


Aie (n; — mi)! 


| Qi e rou - n + (ni m)(z, - n, i212. 
j=1 
MUTER I AII I FUR ZCRBUSCOHPLARIRSEDR L — LiLo 
2 
L(e1,02) = [ [ hie 9:79? /oi2m. (6) 
i-1 


联合 对 数 似 然 函 数 为 


2 2 2 
InL(o1,02) = una 一 X 2m; Inc; 一 Xx Qi/202, 
i=1 i=1 i=1 


olnL 2m1 Qi 
de, 5 to i 
L 


第 1 期 胡 俊 梅 等 ， 双 参数 Rayleigh 分 布 环境 因子 的 经 验 Bayes 估计 及 应 用 


87 


令 
olnL - 0 OinL | 


Oo1 ; Oo» 
从 (7), (8) 式 解 出 o1, o2 的 极 大 似 然 估计 GEX Gir, 02) 分 别 为 


iur = (Qi/2mi)?, G2mr = (Q2/2mm2)12. 


由 此 得 的 极 大 似 然 估计 (i829 kmr) 为 


kumur = 2MrL /1ML = (Qomi/Qimz) 2. 


4 kh 已 知 时 ， 环 境 因 子 k 的 经验 Bayes (EB) 估计 
由 文献 [ 9， 可 取 o; 的 先 验 密 度 为 
Tai) = bo; je 720?, i= 1,2, 
其 中 6b > 0 为 超 参 数 。 由 (5), (10) RE c; 的 后 验 密度 为 


(Qi Fi) aus SET 


rla 19) = 3mm, t 7* 


由 k = 02/01 得 大 的 后 验 密度 为 


十 co 
fk)- J o172(ko1 mi (o1)doy = Ck? [ba + Qo + K* (0, Q)] 017919. 
0 


其 中 
C = 2(Qi 4- b)" *1(Qs + b3)??*1 /B(m +1, mı + 1). 


B(z,y) 为 Beta 函数 。 
在 平方 损失 下 ， 上 大 的 Bayes fiit ke 为 


nos / Da 
0 


il 


上 式 不 能 求 出 精确 解 ， 可 用 数值 积分 求 出 。 
下 面 用 极 大 似 然 方法 估计 超 参 数 b;, i = 1,2 


Too l _ 2bi(zi — u) 
fea - [ A ur am 


十 co bj 
1 — F(z;) = I f(zi)dz; = (zi ug Fb 


似 然 函 数 (4) 式 变 为 


Li(zi | bi) = CR Ius ) | s (asi 


Too 
C Í k?™ +? [ba + Q2 + k?(by + Q1)] imita 
0 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


88 Roo 


学 第 27 卷 


将 (13), (14) 式 代入 (15) 式 得 


Li(zi | bi) = Di II [mE x 
j21 t7 


e uy + bp 


其 中 D; = [(n; 7 m;)!] 127*n! ID za - p). 


bi 


pm 
(Zim, a uy? 十 bi 


| 


? 


mi b; 
In L;(zx; | bi) = ln Di ilnb; — 2 I ij — A b, i — mi) Iln ————————. 
n Li(z; | bi) = ln Dj + mi ln 2i [s LY t bi] + (n mol c EN 
由 
1n Li(xi]b;) z Mi - Mi 2 Ni — Mi T Ni — Mi E 
ðb; j3 bi pen (zij — u)? + bi bi (Zim, — u)? 4- bi ERE 
得 
mm 
bi (cy) +bi (Eim; — 4)? + bi 
令 
Mia 2 n;—m; 
hoje ou) oo o duc Hs e 
1(b:) bi 2(bi) 2 G cu rh Cm ah 
则 有 
hi(bi) > 0， hib) = —3 <0, hin(b:;) = > 0. 


2 


T h: (bi) 在 (0, +oo) 上 为 严格 单调 递减 凹 函 数 


ha(bi) > 0, h3l(b;) 一 一 
2-1; 


leg uy BB [ze ix 


4 


hat(b,) 一 和 = [cw E ME + b;] 


a ha(bi) 在 (0, 十 co) 上 为 严格 单调 递减 止 函 数 。 


im ha (b:) = 
b—oo ha(b) 


故 方程 D 


2 


Ni 


bi Z (ty =u)? + 


具有 唯一 解 ， 且 和 迭代 公式 为 


pot? 一 二 ie 


2 Ni — Mi 


0, 


2(ni = mi) 


——————s———»0.. 
[((Zim, — u)? + b]? 


X 


Ni 
Ni + Mi 


<1 


, 


十 Ni 一 Mi 
bi (Zim; — H)? + bi 


, p-112,-., (16) 


D 2 +m 
(xi, —u)2 zb 
Ji (Cup) C (xim, - i)? eti 


id (16) 3E JE EGARTCEDS bi b; I (12) 式 即 可 求 出 及 的 EB 估计 ken. 
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5 ”随机 模拟 


考虑 应 用 环境 因子 将 较 良好 环境 下 的 试验 信息 折合 到 较 严 酪 环境 下 ， 则 环境 2 下 的 第 ;个 单 
元 的 信息 转化 为 环境 1 的 信息 是 x2;/k， 再 根据 定数 截 尾 的 特点 将 这 m1 + ma 个 单元 信息 进行 
排序 得 到 zi < zz <- < Em, m= M + Mmo 
服从 Rayleigh 分 布 的 部 件 的 可 靠 度 函数 和 失效 率 分 别 为 
Raae m dt 17) 


c? 


截 尾数 据 下 Rayleigh 分 布 参数 的 极 大 似 然 估计 为 


m 1/2 
2; — 21)? + (n = m)(zs — 21)? 


ÉML-—€1, OML = ;N= n +n (18) 


2m 


YE âme, ôme 代入 (17) 式 可 得 到 考虑 环境 因子 的 可 靠 度 和 失效 率 估计 。 

以 上 已 经 得 到 了 Rayleigh 分 布 环境 因子 的 极 大 似 然 估计 和 EB 估计 ， 现 用 Monte-Carlo 方 法 
进行 模拟 比较 。 具 体 步 又 如 下 。 

1) 产生 ni 个 服从 U(0,1) 分 布 相互 独立 的 随机 样本 六 , Yo, , Yn, i 二 1,2; 

2) £3Eu—10,01—4,02—8, $ 


1/2 


zij = |- 20?ln(1 — Y;)) +u, j=1,2,.. ,nm 


则 Til, Ti2, ini 是 服从 Rayleigh 分 布 F(2zi;ci) 的 Tli 个 样本 ; 

3) ”预先 给 定 截 尾数 据 m;， 产 生 样本 zi = (mas Liz, , Limi) H (9) 式 即 可 求 出 ku: 

4) 利用 (16) ok H5, Eb; RA (12) 式 即 可 求 出 RpBp， 由 (17), (18) 式 得 到 在 严酷 环境 下 
考虑 环境 因子 的 用 ,56, Ê, 从 和 不 考虑 环境 因子 的 瑟 o, R, F。 以 上 步骤 随机 模拟 1000 次 ， 计 算 
出 不 同样 本 mw MRE m 下 ， 极 大 似 然 估计 和 EB 估计 的 均值 和 均 方 误差 (MSE)。 环 境 因子 的 
真 值 k = 2。 模 拟 结果 列 于 表 1、 表 2、 表 3 中 。 


表 1: 天 的 MLE 估 计 和 EB 估 计 的 均值 和 均 方 误差 


(n1, n2) (mi, ma) kmt kes MSE(kmz) MSE(kgBp) 
(10,10) (7,7) 2.2146 2.0235 0.5329 0.0148 
(15,12) (11,10) 2.1248 2.0075 0.4636 0.0048 
(20,18) (16,13) 2.1756 1.9904 0.4075 5.5876e-4 
(30,25) (20,15) 1.9507 1.9937 0.2971 9.4102e-4 


(35,30) (28,25) 1.9915 2.1020 0.1173 0.0095 
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A2: 严酷 环境 下 p, al 估计 的 均值 和 均 方 误差 
(mn2) (mi, m2) A(MSE) A(MSE) 51(MSE) ôı (MSE) 
(10,10) (7,7) 12.1251(1.5762) 11.0716(0.7489) 2.7136(1.3756) 3.3561(0.7325) 
(15,12) (11,10) 13.1014(1.7492) 10.5847(0.6543) 2.7531(1.4869) 3.3174(0.6500) 
(20,18) (16,13) 13.0910(2.3549) 10.1602(0.1817) 2.7494(1.5593) 3.5800(0.4213) 
(30,25) (20,15) 11.0955(1.0082) 11.0544(0.8629) 3.0054(1.0052) 3.4250(0.2105) 
(35,30) (28,25) 11.0915(0.9811) 10.3256(0.3858) 2.9851(1.3247) 3.6877(0.6692) 
表 3: 严酷 环境 下 R, > 估计 的 均值 和 均 方 误差 (z — 15, R(15) = 0.4578, r(15) = 0.3125) 
(mna) (mi,m2) R(MSE) R(MSE) F(MSE) f(MSE) 
(10,10) (7,7) 0.6259(0.2682) ^ 0.3914(0.0067) ^ 0.6254(0.3429)  0.0030(0.0526) 
(15,12) (11,10) 0.6524(0.4253) ^ 0.5443(0.0172) ^ 0.5253(0.2285) 0.4000(0.0127) 
(20,18) (16,13) 0.7251(0.5387) ^ 0.5524(0.0118)  0.4987(0.1900) ^ 0.2580(0.1155) 
(30,25) (20,15) 0.6211(0.4201) 0.4784(0.0905) ^ 0.3127(0.0623) ^ 0.2014(0.2259) 
(35,30) (28,25) 0.5449(0.1175)  0.5763(0.1349) ^ 0.1255(0.3826) 0.3157(0.0673) 


从 表 1 可 以 看 出 ， 经 验 Bayes 估计 的 MSE 比 极 大 似 然 估计 的 MSE 小 ， 并 且 样 本 容量 越 小 ， 
截 尾 程度 越 高 ， 效 果 越 明显 。 从 表 2 和 表 3 可 以 看 出 ， 考 虑 环境 因子 的 参数 估计 和 可 靠 性 指标 
估计 要 优 于 没有 考虑 环境 因子 时 的 参数 估计 和 可 靠 性 指标 估计 。 


6 ”结论 


在 可 靠 性 评定 中 ， 环 境 因 子 作为 一 种 扩大 样本 量 的 途径 ， 对 提高 产品 的 可 靠 性 评估 结果 有 
重要 的 意义 。 本 文 利用 经 验 Bayes 法 ， 给 出 了 Rayleigh 分 布 环境 因子 的 经 验 Bayes 估计 ， 遂 过 
理论 分 析 和 模拟 结果 比较 ， 得 出 了 Rayleigh 分 布 环境 因子 的 经 验 Bayes 估计 优 于 极 大 似 然 估 
计 。 现 有 的 环境 因子 的 估计 方法 ， 如 极 大 似 然 估 计 等 ， 是 在 大 样本 理论 的 基础 上 提出 的 。 而 本 
文 提出 的 经 验 Bayes 法 可 用 于 小 样本 试验 的 统计 分 析 ， 并 且 可 充分 利用 现场 试验 数据 ， 因 此 本 
文 提 出 的 方法 更 符合 工程 实际 。 
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Empirical Bayes Estimation and Application for the Environmental 
Factor of Two-parameter Rayleigh Distribution 


HU Jun-mei, SHI Yi-min, GAO Ni, QIN Xiao-qiong 


(Department of Applied Mathematics, Northwestern Polytechnical University, Xi'an 710072) 


Abstract: When the location parameter is known, we discuss the maximum likelihood estimation 
(MLE) and empirical Bayes estimation (EBE) of environmental factor under a two-parameter Rayleigh 
distribution. Furthermore, we estimate the reliability performance by using the estimation of envi- 
ronmental factor. Finally, based on the Monte-Carlo simulation, the simulating results show that the 
accuracy of the EBE is better than that of the MLE and that the reliability assessment is better by 
considering the environmental factor. 

Keywords: Rayleigh distribution; environmental factor; empirical Bayes estimation; reliability assess- 
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